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Abstract 

The andante regime of scalar field dynamics in the chaotic inflationary Universe is 
defined as the epoch when the field is rolling moderately slowly down its interaction 
potential, but at such a rate that first-order corrections to the slow-roll approximation 
become important. These conditions should apply towards the end of inflation as the 
field approaches the global minimum of the potential. Solutions to the Einstein-scalar 
field equations for the class of power law potentials 1^(0) oc 0^" are found in this regime 
in terms of the inverse error function. 
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1 Introduction 



A study of the evolution of self-interacting scalar fields in the early Universe is important 
for a number of reasons. Firstly, in addition to offering a possible resolution to some of 
the fundamental problems of the hot big bang model, the inflationary scenario provides a 
causal mechanism for generating adiabatic density perturbations These may produce 
anisotropics in the cosmic microwave background and act as a seed for galaxy formation 
via gravitational instability. During inflation the Universe is dominated by the potential 
energy V(4>) associated with the self-interactions of a quantum scalar field (p. If the field is 
initially displaced from the global minimum of the potential and the potential is sufficiently 
flat, the scalar field will evolve very slowly towards the true vacuum state. The potential 
therefore behaves as an effective cosmological constant and introduces a negative pressure 
into the Universe that drives the accelerated expansion. Secondly, many natural extensions 
to General Relativity may be expressed in terms of Einstein gravity minimally coupled 
to a self-interacting scalar field after a suitable conformal transformation on the metric 
tensor. Two classes of theory that have received much attention in recent years are scalar- 
tensor theories of gravity and higher-order theories, where the gravitational lagrangian is 
an analytic function of the Ricci scalar R. 

However, only a limited number of exact solutions to the isotropic Einstein-scalar field 
equations have been found to date. These have recently been classified by Barrow in terms 
of the potential V = Vocf)'^ ex.p{A(f>^), where {Vq, A, M, N} are constants Special cases 
include V = constant corresponding to exponential expansion, power-law inflation from 
an exponential potential {N = 0, M = 1) [§, |5[ |^, 0] and intermediate inflation from a 
combination of power-law potentials {N < 0,M = 0) 0, |8|. Solutions corresponding to 



potentials leading to hybrid inflation have been presented in Refs. ||lO| , 11] and exact 
solutions have also been found for hyperbolic and trigonometric potentials p^ . 

Consequently, it is common practice to invoke the slow-roll approximation. This assumes 
that the kinetic energy of the scalar field is negligible relative to its potential energy and 
that the dominant term in the scalar field equation is due to friction arising from the 
expansion of the Universe. The system is therefore reduced to a set of coupled, first-order 
differential equations. The slow-roll approximation is usually valid during the initial stages 
of the infiationary expansion, but as the field rolls towards the minimum, the approximation 
inevitably breaks down at some point. Moreover, inflation arises whenever the strong energy 
condition is violated and this does not necessarily require the slow-roll approximation to be 
valid. 

Therefore, it is important to develop analytical techniques that allow deviations from the 
slow-roll regime to be accounted for and in this paper we discuss how this may be achieved. 
After summarizing the general features of scalar field dynamics in Sec. II, we illustrate how 
corrections to the slow-roll approximation may be included in Sec. Ill and we derive the 
corresponding field equations. In Sec. IV these equations are solved in parametric form 
for the class of power law potentials V{(j)) oc (j)"^^ in terms of the inverse error function. 
The rapid increase in the scale factor during inflation is naturally explained by using the 
properties of this function. 



1 



2 Scalar Field Dynamics 



Considerable progress in determining scalar field dynamics has been made recently by treat- 
ing the scalar field as the dynamical variable, since this reduces the field equations to a set 
of first-order, non-linear differential equations 0, |5|, ^ U, |l5|. The four-dimensional 
action for Einstein gravity minimally coupled to a self-interacting scalar field is 



S 



-a 



R 

2^ 



(1) 



Svrmp^ and mp is the Planck mass. We choose units such that 



where g 

c = h = 1. If the Universe is spatially isotropic, closed and flat with a world-interval 
ds'^ = -N'^{t)df + e^"(*)[d3;2 + dy"^ + dz^] and lapse function N{t), the Arnowitt-Deser- 
Misner (ADM) action is 



S 



dtUNe 



3a 



3 1 T..^^' 



(2) 



where U = J d^x. is the comoving volume of the Universe and a dot denotes differentiation 
with respect to cosmic time t [16|. The classical Hamiltonian constraint Ti = implies that 



(3) 



the action satisfies the Hamilton-Jacobi equation 



K 



+ 



dS\ 



j + 2C/2e6"y((/>) = 0. 



.da J 

The classical dynamics of this model is determined by the real, separable solution 

5= -^C/e3"i/(</>), 
where H{(j)) satisfies the differential equation |0, [s], 



--^V(<^). 



(4) 



(5) 



This equation is equivalent to the 00-component of the Einstein field equations and therefore 
represents the Friedmann equation. In the conformal gauge = 1 the momenta conjugate 
to a and are pa = dS/da = — 6K~^C/e^°d and = dS/dcf) = Ue^"(j), respectively. 
Substitution of ansatz (H) into these expressions implies that 



a. 



2 dH 

d(j) 



(6) 



and it follows that H{(j)) is the Hubble expansion parameter expressed as a function of the 
scalar field cp. The general solution to the field equations (^) and @ may be expressed in 
terms of quadrature with respect to this function: 



a 



K 

~2 



b'H{<p',p) 



fdH{(P',p) 



(7) 
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where {ai,(j)i) are arbitrary constants and p is a parameter associated with each solution 
p) ofEq. (ID I, I- 

In principle the general path of the Universe in the minisuperspace (a, (j)) is uniquely 
determined once the functional form of H{(f>) is specified and this suggests that it should be 
be viewed as the fundamental quantity in the analysis. Unfortunately, however, it is rather 
difficult to find exact solutions to Eq. (^) that go beyond the slow-roll approximation. In 
view of this it has been suggested that one could begin the analysis by specifying H{(p) 
1 11, 17]. An alternative approach is to generate exact solutions by viewing the expansion 
parameter as the effective time coordinate |[T^ . The drawback with these approaches, 
however, is that the exact solutions do not necessarily correspond to realistic potentials. 
Moreover, since it is the potential, and not H{cj}), that is specified by the particle physics of 
the model, one should aim to solve the model by specifying the potential. In the next section 
we show how an approximation to Eq. (^) can be made when the slow-roll approximation 
is not valid. 



3 Beyond Slow-Roll: The Andante Regime 

It proves convenient to introduce the rescaled quantities |Q] 

,2 \ 1/2 ^3\V2 



U 



-V 



X 



since Eq. (P) then simplifies to 



(8) 



(9) 



Deviations from the slow-roll regime have been analyzed previously by Salopek and Bond |5[. 
In order to derive a more accurate solution they substituted the zeroth-order approximation 
HiQ\ = u into the right-hand side of Eq. (0) : 



(1) 



(0) 



1 + 



dx 



(10) 



and further accuracy is achieved by including higher-order terms. 

However, deviations from the slow-roll approximation may be studied more quantita- 
tively in terms of the dimensionless parameters 



3x 



i.2 



x^ + Qn^ 



and 



H± 



H" 
'IT 



1 /3V/2 , 



e - - - 



H'" 



3\i/2 

r] , 



(11) 



(12) 



where a prime denotes differentiation with respect to the dimensionless field x For 
convenience we consider i; > in this work, which implies that ^/e = —\/^H'/H. Modulo 
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a constant of proportionality, e is a measure of the field's kinetic energy relative to its total 
energy, whilst r] measures the field's acceleration relative to the amount of friction acting 
on it due to the expansion of the Universe. We may therefore refer to these quantities as 
the energy and friction parameters, respectively The slow-roll approximation to scalar 
field dynamics corresponds to ^ 1 iii this regime the energy and friction 

parameters reduce to the slow-roll parameters introduced in Ref. |20|. It is straightforward 
to show that inflation occurs if e < 1 and the end of inflation can be defined precisely by 
the condition e = 1 |, |, |l7|. 

Eq. dll) can be written in a very illuminating form by introducing the parameter 



e 



1. 



It follows that 

1 _ e _ u"^ 
^ ~ 3 ~ ~lP 
and differentiation of this equation with respect to x implies that 



u' _ 1 v' 
u V v{v'^ — 1) 



(13) 
(14) 

(15) 



A comparison of Eqs. (13) and (^) implies that the 'steepness' of the potential, as defined 



by Turner [21|, may now be expressed directly in terms of the energy and friction parameters: 



u 

u 



1 - e/3 



(16) 



The rate of change of steepness in the potential can also be written with these parameters 
and has the form 



1 d^V 



u I u 



r?[l - 77/3] + e[l - C/3] 
3[1 - e/3] 



(17) 



Hence, the first slow-roll condition corresponds to u'/u ~ — ^e/3 <C 1, i.e. ii(x) « u(x). 

We may consider first-order departures from the slow-roll regime by expanding the terms 
in the square bracket of Eq. (|l^) to first-order in e and 77. With the help of Eq. ( p!^ ) we 
find 

v! 1 v' ( v' 



(18) 

U V ' 

and the dymanics of the scalar field in this regime of parameter space is then determined 
by this equation. Eq. ( [T8| ) may also be derived by expanding the right-hand side of Eq. 
( [TsD as a geometric progression for v > 1: 



u 
u 



1 V 

-- + -. 

V v 



1 1 
l + ^ + ^ + .-- 



(19) 



It is seen, therefore, that accounting for first-order deviations from the slow-roll regime is 
equivalent to including the first term in this series. 
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For a specific choice of potential the solution to Eq. ([l8|) will determine v{x) and 



the functional form for H{x) then follows directly from Eq. (14). Eq. (18) describes 
the dynamics of a scalar field that is evolving moderately slowly down its self-interaction 
potential but at a sufficiently fast rate that the usual slow-roll approximation is no longer 
valid. Consequently we refer to this region as the moderately slow, or andante, regime 
of scalar field dynamics^. The slow-roll limit corresponds to the regime where the energy 
and friction parameters are negligible with respect to unity and the andante regime applies 
when these parameters are small but finite. Expressions for the amplitudes of the scalar 
and tensor fluctuations when these latter conditions apply have been derived previously by 



Stewart and Lyth |22|. 



An alternative form of Eq. (18) may be derived by viewing the volume of the Universe 
as the effective dynamical variable. Since the volume is a monotonically increasing function, 
it is a suitable choice for the 'time' coordinate. It follows from the general solution (0) that 
the number of e-foldings of expansion that occur as the scalar field rolls from some initial 
value Xi to a value x is given by 

Ina oc N{x) = t dt'H{t') = -- ^ dxH{x') ( ^EMl) 
Jti 3 Jxi \ dx' ) 

= - r dx'v{x'). (20) 
3 

Hence, the ratio zj Zi = (a/ai)^ = exp(3A^) represents the fractional increase in the volume 
of the universe that occurs when the field rolls from Xi to x and it follows immediately from 
the definition (|^) that 

- = -• (21) 

dz zv 

If we now define a new function m{z): 

m{z) = zx{z), (22) 
it can be shown, after differentiating with respect to z and substituting in Eq. (^), that 

dm 1 

— = X + - (23) 
dz V 

and 

d'^m u' 

m = —X — . (24) 
dz'^ u 



We shall show in the following section that Eq. (24) is exactly solvable when the potential 
is given by V{4>) oc 0^'^. 

^Similarly the slow-roll regime could be referred to as the 'largo' regime, whilst the epoch after inflation 
that arises when the field undergoes rapid oscillations about the minimum would represent the 'allegro' 
regime. 
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4 Power Law Potentials 



In this section we consider the class of power law potentials 



T/(0) = 



n = constant. 



(25) 



where A is the coupling constant. In the chaotic inflationary scenario, the andante regime 
should apply during the final stages of the inflationary expansion when the scalar field 
approaches the global minimum of its potential. In this region it is useful to consider the 
general inflationary potential as a truncated Taylor series expanded about this minimum. 
It is therefore a good approximation to assume the potential has the power law form ( p5|) 
with positive-definite n. However, the class of potentials with n < is also interesting and 
potentials of this form may be important if the Universe contains a nonvanishing vacuum 



energy at the present epoch |23]. Moreover, they arise in generalized scalar-tensor theories 
when the Brans-Dicke parameter is viewed as a truncated Taylor series in the dilaton field 

When the potential is a power law, Eq. (|2^ ) admits the first integral 



dm 
dz 



lb (c — 2?7-lnm) 



1/2 



(26) 



where c is an arbitrary constant. For consistency we require m > in Eq. (|26|), and since 
we assumed x > 0, this implies that both x and z must be negative if n > and they must 
be positive if n < 0. However, \z\ provides a measure of the amount of inflation that occurs 
and it varies much more rapidly than the scalar field x. Thus, the evolution of the function 
m is dominated by z, which implies that m increases as \z\ increases. Consequently, the 
negative square root should be chosen when n > and the positive square root corresponds 
n < . 

Eq. (|26| ) is solved exactly in terms of the error function: 



z = b + 7erf 



T 1 - - Hzx) 



where b is the second integration constant and 
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2nJ 



(27) 



(28) 



This solution may be derived from the identity d[erf{y)]/dy = (2/y^) exp(— y^) |25] 
Eq. ( |27| ) can be inverted to yield 



1 



X{z) =e^/(2n)_g^p 

z 



lerf 



z-b 

7 



(29) 



where lerf (y) is the inverse error function and the expression for v{z) follows by substituting 
Eqs. (H) and (H) into Eq. (|2|): 



— = —X — \/2nIerf ( - — - 
V V 7 



(30) 
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Finally we find z{x) by substituting Eq. (30) into Eq. (|2 



z{x) 



,c/(2n)l 



exp 



1 / ^ 1 

2n \ v{x) 



(31) 



These solutions exhibit the correct qualitative behaviour that one would expect for a 
scalar field rolling down a polynomial potential. Figures la and lb illustrate the behaviour 
of the functions x{z) and v~^{z) for n = 1, c = and b = —7 and similar results are found 
for others choices of these constants. A relatively small change in the value of the scalar 
field in the range —30 < x < —10 results in a huge change in the value of z by a factor 
~ lO^^''. The origin of this behaviour is traced to the properties of the inverse error function. 
Ierf(?/) is an odd function of y and is undefined for \y\ > 1. For positive arguments it is a 
monotonically increasing function with lerf (0) = and lerf (1) = 00. However, the slightest 
deviation of the argument from unity results in a sharp decrease in the value of the inverse 
error function. For example, Ierf(l — 10"^'^"'^) = 47.9 and Ierf(l — lO^^'^'') = 15.1. This 
implies that lerf(y) is a relatively slowly varying function in the range < y < 1 — 10^^'^^ 
and consequently z may change by a factor lO^'' without there being a significant change 
in the value of the scalar field. This feature leads to a rapid growth in the scale factor for 
a very small change in the value of the scalar field. 

In the limit ^ the functional form of the solution ( pl| ) implies that the scale 
factor may be expressed directly in terms of the scalar field as 



From this expression it follows that Eq. (20) may be employed to evaluate H{x): 



n 

X'^ 



(32) 



(33) 



and if one expands the right-hand side to first-order, the result is equivalent to Eq. (10) 



Hence, the Salopek-Bond approximation Q is recovered from the parametric solution (29)- 
(|30|) in the limit that v diverges. 

Similar arguments apply when n < 0. Considering negative values of n is equivalent 
to Wick rotating 7 to the imaginary axis, i.e. 7 i^. The error function with imaginary 
argument is related to the imaginary error function, QT:i{iy) = zerfi(y), and it follows that 
solutions (p9|) and (|30|) take the form 



x{z) 



e'=/(2'^)iexp 



lerfi 




and 



2|n|Ierfi 




(34) 



(35) 



respectively. Figures 2a and 2b illustrate the evolution of x{z) and v^^{z). The De Sitter 
solution is the attractor at infinity for these models and the intermediate infiationary so- 
lution is recovered at large x [^. Solutions ( |34|) and ( |35|) illustrate analytically how these 
asymptotic solutions are reached. 
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5 Conclusion 



In this paper we have considered deviations from the slow-roh approximation by including 
the first-order contributions from the kinetic energy and acceleration of the field. This allows 
a more accurate analysis of scalar field dynamics to be performed. Such improvements are 
expected to be relevant towards the final stages of inflation when the scalar field is close 
to the global minimum of its potential. In view of this we considered a class of power 
law potentials and an exact parametric solution to the field equations was found in terms 
of the inverse error function. This function exhibits some very interesting properties and 
it was shown how the rapid growth in the scale factor of the Universe during inflation is 
naturally explained by employing the properties of this function. Since a large increase in 
the scale factor for a very small change in the value of the scalar field is a generic feature 
of the inflationary scenario, we conjecture that the inverse error function may arise in the 
general solution to the field equations. This possibility is currently under investigation. The 
solutions presented here should improve our understanding of how the Universe moves out 
of the inflationary epoch and into the reheating phase. It will be interesting to investigate 
whether solutions to Eqs. ( [l^ ) and (p^) can be found for other potentials. 



Note: After this work was completed we received a preprint by Liddle et al. that also 



discusses the slow-roll approximation though with a different method of analysis [26|. 
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Figure Captions 



Figure 1: (a) A plot of the solution (29) representing the evolution of the scalar field x 
with respect to — log^gz for the quadratic potential u{x) oc x (n = 1). The numerical values 
of the integration constants are specified to be c = and h = — (7r/2)^/^. (b) The solution 
for the same choice of parameters as in Figure la. 



Figure 2: (a) Illustrating the solution (^) for the decaying power law potential u cc 
{n = —1), where b = c = and 7 = (7r/2)^/^. (b) The corresponding solution (^). The 
solutions approach the intermediate inflationary solution and the De Sitter solution is the 
attractor at infinity. Hence, there is no exit from inflation for these potentials. 
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